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Abstract. The "dark energy" problem is investigated in the framework 
\ of the Poincare gauge theory of gravity in 4-dimensional Riemann-Cartan 

\ space-time. By using general expression for gravitational Lagrangian ho- 

mogeneous isotropic cosmological models with pseudoscalar torsion func- 
tion are built and investigated. It is shown that by certain restrictions on 
indefinite parameters of gravitational Lagrangian the cosmological equa- 
tions at asymptotics contain effective cosmological constant and can lead 
\Q [ to observable acceleration of cosmological expansion. This effect has geo- 

metrical nature and is connected with space-time torsion. 



1 Introduction 



Cosmological researches fulfilled during several last decades were fruitful for theoretical and 
observational cosmology. The creation of inflationary paradigm, which permits to solve a 
number of problems of standard cosmological scenario, in particular, to explain the homo- 
geneity and isotropy of the Universe at initial stages of cosmological expansion, is important 
achievement of the early Universe cosmology [1]. The discovery of the acceleration of cos- 
mological expansion at present epoch is the most principal achievement of observational 
cosmology [2]. By using Friedmann cosmological equations in order to explain accelerating 
cosmological expansion, the notion of dark energy (or quintessence) was introduced in cos- 
mology. According to obtained estimations, approximately 70% energy in our Universe is 
connected with some hypothetical form of gravitating matter - - "dark energy" - - of un- 
known nature. Recently many investigations devoted to dark energy problem were carried 
out (see review [3]). According to widely known opinion, the dark energy is associated with 
cosmological term. If the cosmological term is connected with the vacuum energy density, it 
is necessary to explain, why it has the value close to critical energy density at present epoch. 

The present paper is devoted to investigation of the "dark energy" problem in the frame- 
work of the Poincare gauge theory of gravity (PGTG), which is natural generalization of 
Einsteinian general relativity theory (GR) by applying the gauge approach to theory of grav- 
itational interaction. According to PGTG the physical space-time possesses the structure of 
Riemann-Cartan continuum with curvature and torsion. By using gravitational Lagrangian 
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of PGTG in general form including both a scalar curvature and various invariants quadratic 
in the curvature and torsion tensors (see below), isotropic cosmology was built and investi- 
gated in a number of papers (see [4-6] and references herein). As it was shown, the PGTG 
permits to solve the problem of cosmological singularity of GR: all solutions of generalized 
cosmological Friedmann equations (GCFE) for homogeneous isotropic cosmological models 
(HICM) deduced in the framework of PGTG are regular in metrics, Hubble parameter , 
its time derivative, if gravitating matter satisfies certain physically acceptable condition at 
extremely high energy densities and pressures [4]. It is a consequence of gravitational re- 
pulsion effect, which takes place at extreme conditions and is connected with space-time 
torsion essentially [6]. Properties of discussed HICM in PGTG coincide practically with 
that of GR at sufficiently small energy densities, which are much less in comparison with 
limiting (maximum) energy density for such models. By including cosmological term of cor- 
responding value to the GCFE, we can obtain regular cosmological solutions with observable 
accelerating expansion stage. However, like GR, the problem of dark energy in such theory 
is not solved. 

From geometrical point of view, the structure of HICM in PGTG can be more complicated 
in comparison with models describing by GCFE. Really, in the case of homogeneous isotropic 
models the torsion tensor S x fiv = —S x Ufl can have the following non-vanishing components 

[7, 8]: S\ = S 2 20 = S% = ^123 = S 231 = S 312 = S 2 (t) Rr sine, where S 1 and 

S 2 are two functions of time, spatial spherical coordinates are used. The functions Si and S 2 
have different properties with respect to transformations of spatial inversions, namely, the 
function S 2 (t) has pseudoscalar character. The GCFE follow from gravitational equations 
of PGTG for HICM together with S 2 = 0. Obtained physical consequences of GCFE have 
principal character. However, it is necessary to note that gravitational equations of PGTG 
for HICM have also other solution with non- vanishing function S 2 . The HICM with two 
torsion functions are studied below in the frame of PGTG in connection with the dark energy 
problem. In Section 2 gravitational equations of PGTG for HICM with two torsion functions 
are obtained and investigated. In Section 3 the asymptotics of cosmological solutions for such 
models is analyzed. 

2 Homogeneous isotropic models with two torsion func- 
tions in PGTG 

At fist let us mention some general relations of the PGTG. Gravitational field is described 
in the frame of PGTG by means of the orthonormalized tetrad h 1 ^ and anholonomic Lorentz 
connection A lk ^ (tetrad and holonomic indices are denoted by latin and greek respectively); 
corresponding field strengths are torsion S\ v and curvature F tk ^ tensors defined as 

F;t = 2^ + 2^4 M . 
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We will consider the PGTG based on the following general form of gravitational Lagrangian 



C G = h [f F + F a ^ (/i F a ^ + f 2 F a ^ v + h F^p) + F^ (/ 4 F, v + f 5 F UfJ ) + / 6 F 2 

+5^ ( fll + a 2 S^) + a 3 5 W' 3 ] , (1) 

where h = det F^ = F a ^ au , F = (i = 1, 2, . . . , 6), a k (k = 1, 2, 3) are indefinite 

parameters, fo = (16irG)~ 1 , G is Newton's gravitational constant. Gravitational equations 
of PGTG obtained from the action integral I = f {C g + C m ) d^x, where C rn is the Lagrangian 
of matter, contain the system of 16+24 equations corresponding to gravitational variables 
^VandA%. 

Any homogeneous isotropic gravitating system in PGTG is characterized in general case 
by three functions of time: the scale factor of Robertson- Walker metrics R and two torsion 
functions S\ and S 2 . Below by using spherical coordinate system, the tetrad is taken in 
diagonal form. Then the curvature tensor has the following non- vanishing tetrad components 
denoted by means of the sign ' : 

-F23 23 = A 2 , 



with 



b 6i - 

TP 23 

^oi - 


r 02 
r 02 


- TP.K 

— r 03 

~~ ^03 


= A 3 , 


TP 12 
b 12 - 

p..6i - 
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A 2 = 




- 2S,) 2 - 


o2 
D 2> 



F21 = A 4 , 



A 3 = 2(H-2S 1 )S 2 , 

A4 = S 2 + -^5*2, 



(2) 



where H = R/R is the Hubble parameter and a dot denotes the differentiation with respect 
to time. 

The Bianchi identities in this case are reduced to two following relations: 

A 2 + 2H (A 2 - A x ) + + 2S 2 A A = 0, 

A 3 + 2H {A 3 - A 4 ) + 45^4 - 2S 2 A± = 0. 

By using the gravitational Lagrangian (1) and the Bianchi identities (3), we obtain the 
following system of gravitational equations for HICM 

a (H - 50 5 X - 2bS 2 2 - 2f A 2 + 4/ (A\ - A 2 2 ) + 2q 2 (A\ - A 2 A ) = -|, (4) 

a (S 1 ! + 2HS 1 - 5 2 ) - 2bS 2 2 - 2f (2A, + A 2 ) - Af (A\ - A 2 2 ) - 2q 2 (A 2 3 - A 2 ,) = p, (5) 

f[(Ai + i 2 ) + AS, (A, + A 2 )] + g 2 S 2 A 3 + (2/ - Ql ) S 2 A 4 + (/ + |) ^ = 0, (6) 

(i 3 + A^j + AS, (A 3 + AO] - 4/5 2 A 2 - 2 ( 9l + g 2 ) - (/ - 6) 5 2 = 0, (7) 
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where 



2<2i + a 2 + 3a 3 , 

h 



a 2 - ai, 



/ = A + f + /3 + /4 + / 5 + 3/ 6 , 

5i = / 2 - 2/ 3 + / 4 + /s + 6/ 6 , g 2 = 2/x - / 2 , 

p is the energy density, p is the pressure and the average of spin distribution of gravitating 
matter is supposed to be equal to zero. Equations (4)-(5) do not contain high derivatives 
for the scale factor R, if a = 0, moreover, equations (6)-(7) take more symmetric form, if 
2/ = <Zi + Q2- Then the system of gravitational equations for HICM take the following form: 
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- 2b S 2 2 - 2f A 2 + 4/ (A\ - Aj) + 2q 2 (A 2 3 - A 2 4 ) = --p, 

- 2b Si - 2f (2A, + A 2 ) - Af (Aj - A\) - 2q 2 (Aj - Aj) = p, 
f[(Ai + i 2 ) + AS! (A, + A 2 )] + q 2 S 2 (A 3 + A 4 ) + f Q S! = 0, 
A 3 + AA + AS, (A 3 + A,)} - Af S 2 (A ± + A 2 ) - (f -b)S 2 = 0. 



(8) 

(9) 
(10) 

(11) 



The system of equations (8)— (11) together with definition of curvature functions (2) is the 
base of our investigation of HICM below. Note also that the conservation law for spinless 
matter has usual form: 

p + 3H(p + p) = 0. (12) 

From (8)-(9) follows that 



Ax + A 2 



(13) 



By using (13) and the formula following from definition (2) of curvature functions 



A 2 - A 2 



4A 2 5f-4(A 



si ) si 



S 2 + HS 2 



we find from gravitational equations (8)-(9) the following expressions for A\ and A 2 : 



A 1 = 



1 



12/oZ 



p + 3p 
+ 



Ao = 



QfoZ 



a 
~ 2 

3ae/ 
Z 

a 



2\2 



p-3p- 12bSl) 



ae 



[p-3p- 12bSl) S 



HS 2 + S 2 ) ±A(—-Sl)S 2 



p-QbSl + -(p-3p-12bSiy 



3aef 



(HS 2 + S 2 y + A^-Sl ) jS 



(14) 



(15) 
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where Z = 1 + a [(p — 3p) — 12 (b + ef ) Sf], a = g4-, e = — , — = 3a ef . By using the 

/o / to 

conservation law (12), the formula (13) and the following relation obtained from definition 
of A 3 and A 4 



A 3 + A 4 = S 2 + 3HS 2 - ASA 
we obtain from (10) the following expression for the torsion function Si. 



Si 



3a 
AZ 



where 



A(2b-sf )S 2 S 2 -HY 



dp 



(16) 



(17) 



Y = (p + p)(3f p -l) + l2ef Si 



Then by using formulas (13) and (16) we find from (11) the following differential equation 
of second order for the torsion function S 2 : 



S 2 + 3HS 2 + (3H - ASi ) S 2 + AS^ (3H - AS ± ) 



Uo ~ b) 



TT { P -3p-l2bSl)+ , 
6 Jo J 



S 2 = 0. (18) 



The obtained expressions (14)— (15) for curvature functions A 2 and A\ together with their 
definition (2) give the generalization of cosmological Friedmann equations for HICM: 



k/R 2 + (H- 2S 1 ) 2 = —— \p + 6 (f Z - b) S 2 2 + ^( p -3p- 12bS 2 2 ) 2 



GfoZ 



3aef 



(HS 2 + S 2 y + A(^--Sl)s 



, (19) 



H + H 2 -2HS 1 -2S 1 



WoZ 



p + 3p 



+ 



"2^ 
3ae/ 



as 



3p - \2bS 2 ) 



— (p - 3p - 12bS 2 2 ) S 



HS, 



S 2 2 



si 



(20) 



These equations contain the torsion function determined by (17) and the torsion function 
S 2 , satisfying the equation (18). Besides indefinite parameter a determining the scale of 
extremely high energy densities [4], obtained equations contain additionally two indefinite 
parameters: b with dimension of parameter f and the parameter e without dimension. We 
have to analyze all these equations in order to investigate HICM with pseudoscalar torsion 
function in the frame of PGTG. 
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3 Asymptotics of cosmological solutions for HICM with 
pseudoscalar torsion function 



The structure of obtained equations describing HICM with pseudoscalar torsion function is 
essentially more complicated in comparison with the case of HICM with S2 = 0. By using 
the conservation law (12), the cosmological equation (19) can be written in the following 
form: 



H 2 



3 \ 

Z + -^YJ + 3aef S%Z - 6aH 



2\Z+^Y)(2b -cf {l )S 2 S 2 -::f {) S 2 S 2 Z 



+ 36a 2 {2b - efo) S 2 2 S 2 2 + 3aef 



1 

Wo 



k_ 



-si 



s 2 



p-6bSl + ^{p-3p- \2bSiy 



Z + 



k_ 



-s 2 



Z. (21) 



From this equation follows, that cosmological solutions possess extreme points for the scale 
factor R. In fact, if we put that at t — the Hubble parameter vanishes H = 0, then from 
(21) follows the relation for S 2 , its time derivative and energy density, which has to be valid 
at a bounce: 



36a 2 (26-e/o) 2 S 2 20 S 2 20 + 3as f 



4 2 o + 4 
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~ Wo 
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D2 °20 
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Po 



QbS lo + ^ (PO 



3p 



ms 2 20 y 



Z . (22) 



(Values of various quantities at extreme point are denoted by index 0). By using the following 
equation of state p = p and by choosing some initial values for £20, 5*20! Po a t a bounce, 
particular cosmological solution for flat model {k = 0) presented in Fig. 1 - Fig. 2 was 
obtained numerically. We see that the pseudoscalar torsion function S2 and the Hubble 
parameter H have some non- vanishing values at asymptotics. Below we analyze the following 
question: by what restrictions on indefinite parameters b and e asymptotical value for the 
Hubble parameter corresponds to observable accelerating cosmological expansion? 

By taking into account that various parameters of HICM have to be small at asymp- 
totics, we see from (18), that if e 1, the pseudoscalar torsion function has the following 
asymptotics: 

q2 _ fo(fo-b) p-3p 

Note, that at asymptotics one uses for gravitating matter the equation of state for dust 
p = 0; then it follows from (23) that b < fo and we obtain at asymptotics: Z — > {b/fo), 
Si — > 0. The cosmological equations (19)-(20) at asymptotics take the form of cosmological 
Friedmann equations with cosmological constant: 



k zj2 1 

& +H = 6b 



P + 



3 {fo - by 

4/ 



(24) 
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Figure 1: The Hubble parameter H(t) in di- 
mensionless units. 




Figure 2: The pseudoscalar torsion function 
S 2 (t) in dimensionless units. 



H + H = 



12b 



p + 3p 



3 (/o - 6)' 
2/ 



(25) 



From Eqs.(24)-(25) we see, that parameter b has to be very close to / . The value of b leading 
to observable acceleration of cosmological expansion depends on the scale of extremely high 
energy density defined by a -1 . If we take into account that the value of |(/ — b) 2 / f = 
\a~ l {l — b/ f Q ) 2 has to be equal approximately to 0.7p cr (p cr = 6/o-^o> is the value of 
the Hubble parameter at present epoch), then we obtain that b = [1 — (2.8p cr a) 1 ^ 2 ]fo. If 
we suppose that the value of a -1 is greater than the energy density for quark-gluon matter, 
but smaller than the Planckian one, then we can obtain the corresponding estimation for b, 
which is extremely close to / . 



4 Conclusion 

As it was shown, the presence of pseudo-scalar torsion function in HICM built in the frame- 
work of PGTG can lead to effective cosmological constant in asymptotics of cosmological 
solutions and to observable accelerating cosmological expansion. The effect of acceleration 
of cosmological expansion in PGTG has the geometrical nature and is connected with space- 
time torsion. From the point of view of considered theory hypothetical form of gravitating 
matter — dark energy or quintessence — is fiction. The further investigation of HICM with 
pseudoscalar torsion function will be to continued. 
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